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We study the heat transfer between two parallel metallic semi-infinite media with a gap in the 
nanometer-scale range. We show that the near-field radiative heat flux saturates at distances smaller 
than the metal skin depth when using a local dielectric constant and investigate the origin of this 
effect. The effect of non-local corrections is analysed using the Lindhard-Mermin and Boltzmann- 
Mermin models. We find that local and non-local models yield the same heat fluxes for gaps larger 
than 2 nm. Finally, we explain the saturation observed in a recent experiment as a manifestation of 
the skin depth and show that heat is mainly dissipated by eddy currents in metallic bodies. 



I. INTRODUCTION 



Near-field radiative heat transfer has been investi- 




H |2£ |23L |24 |2& |2& I22L \2$). Rytov and co- 
workers |l[ showed how to calculate thermal radi- 
ation by introducing fluctuational electrodynamics. 
This theory is based on the introduction of ran- 
dom current densities due to the thermal random 
motion of charges. Their correlation functions are 
given by the fluctuation-dissipation theorem (FDT). 
Cravalho, Tien and Caren [2j and Olivei [3( were the 
first to address heat transfer in the near field, i.e., 
at distances smaller than the peak wavelength At of 
the thermal radiation spectrum. However, they did 
not consider all evanescent waves. Polder and Van 
Hove [H were the first to take into account all the 
evanescent waves by using the formalism introduced 
by Rytov. They found a huge increase in the heat 
flux between two parallel surfaces when the gap dis- 
tance d becomes smaller than At- Rytov QJJ and 
co-workers pointed out that spatial dispersion could 
play a role for small gaps. Volokitin and Persson [l9| 
showed that spatial dispersion could be responsible 
for an increased heat flux in the nanometer-range by 
using an approximation for the non-local reflection 
coefficients. Loomis and Maris [l4| also investigated 
heat transfer between metallic bodies, showing the 
influence of the electrical resistivity. Recently, Mulet 
et al. [13, showed that the radiative heat transfer 
between dielectrics supporting surface phonon po- 



laritons is dominated by the surface wave contribu- 
tion. As a result, the heat flux is monochromatic in 
this case. 

Several experiments have been reported. Tien's col- 
laborators made the first measurements at cryogenic 
temperatures, when the near field starts at hundreds 
of microns. Kuteladze and Bal'tsevitch [lCjl per- 
formed an analogous experiment. Hargreaves [6| was 
the first one to note (at ambient temperature) an en- 
hanced heat transfer over micrometric distances by 
using two parallel plates of chromium. In the end of 
the 1980's, Xu [Ht [Tl| et al. could not confirm this 
effect with an indium needle in front of silver. Re- 
cently, Kittel et al. [I?], H3| showed a large increase 
in the heat exchange between a scanning probe mi- 
croscope (SPM) metallic tip and a planar surface by 
working in the nanometer range. Surprisingly, they 
also found that the increase of the heat flux levels off 
(saturates) at very small scales (a few 10 nm). This 
is in striking contrast with the 1/d 3 dependence to 
the distance d of the density of states close to the 
surface. It is also in contrast with the power laws dis- 
cussed by Pan et al. [18] . This led Kittel et al. to sug- 
gest that the observed saturation at short distances 
be due to a non-local dielectric constant. Very re- 
cently, Narayanaswamy and Chen [3(| measured an 
enhancement of the heat flux at micron distances, 
using dielectric polar material and a setup similar 
to the one used for measurements of the Casimir 
force [29( • Simultaneously, a number of groups tried 
to use proximity-enhanced heat transfer to increase 
locally the number of electric charge carriers. Di 
Matteo et al. [3l| reported an experimental obser- 
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vation in 2001. A number of theoretical papers also 
present heat flux levels [mSmQSIME 
It has also beenpredicted that metamaterials 
electron doping [25| or adsorbates [3!| may enhance 
the near-field heat transfer. 

Although the enhancement of the flux becomes very 
large at distances on the order of a few nanometers, 
most of the published results use a local model of 
the dielectric constant. It has been pointed out that 
non-local effects should affect significantly the life- 
time of a molecule close to a surface El, H2, |43[ . 
This effect has also been studied in the context of 
the Casimir force 0, Hf| . It appears to be a rela- 
tively minor correction. The experimental findings 
of Kittel et al. has revived the interest for non-local 
effects as the saturation observed at short distance 
is a very significant effect. This paper is devoted to 
the analysis of two questions: i) what is the origin of 
the saturation of the flux in the near field ? ii) what 
are the consequences of non-locality in the context 
of near-field radiative heat transfer ? 

In this paper we focus on the heat flux between two 
parallel semi-infinite metallic substrates. We show 
that for a metal the s-polarized (transverse electric, 
TE) contribution is the leading one in the nanomet- 
ric regime when using local optics. Indeed, the con- 
tribution of the familiar 1/d 2 divergence at short dis- 
tances due to p-polarized waves becomes the lead- 
ing contribution only below 0.1 nm. The satura- 
tion of the s-polarized contribution is similar to the 
experimental behaviour reported by Kittel [27|, so 
that non-local corrections do not seem necessary. 
To further investigate this issue, we compute the 
near-field radiative heat transfer using two non-local 
models: the Lindhard-Mermin model based on the 
random phase approximation and its approximation 
in the Boltzmann-Mermin model. Both longitudinal 
and transverse non-local dielectric constants are in- 
cluded. We find that a local calculation agrees well 
with the non-local ones at gap distances larger than 
2 nm. We finally discuss the physical mechanism re- 
sponsible for the saturation. We show that it is due 
to the magnetic fields that generate eddy currents. 



II. NEAR-FIELD RADIATIVE HEAT FLUX 
USING A LOCAL DIELECTRIC CONSTANT 



We start the section by summarizing the derivation 
of the heat flux between two parallel semi-infinite 
bulks. We do not consider any roughness or tilt 



between the surfaces. Both semi-infinite media are" 
assumed to be in local thermodynamic equilibrium 
with temperature T\ and T 2 . This allows to derive 
the energy radiated by random currents in medium 
1 at temperature T\ and absorbed in medium 2 and 
vice versa. The model can be extended to inhomo- 
geneous temperature profiles provided that the tem- 
perature variation across a distance of the order of 
the skin depth is negligible. The flux per unit area 
is given by the normal component of the Poynting 
vector, 



(1) 



where the position r can be taken at the center of 
the gap z — and (• ■ ■) denotes a statistical average. 
Derivations can be found in many articles [1, [TJ, Q, 
[l6l [H, [H, [5i| and will not be repeated here. 
The final form of the heat flux is 
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where d is the distance between the two interfaces, 
r^ m the reflection factor at the interface between 
medium m and vacuum (medium 3) for a wave with 
wave vector K parallel to the surface and polariza- 
tion a = s, p. The wave vector 



7,„ = ^/e m ui 2 /c 2 - K 2 =i m + i-y'm 
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describes the propagation across medium m, c is the 
speed of light, and 
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is the monochromatic specific intensity of blackbody 
radiation with h, ks the Planck and Boltzmann con- 
stants. We now discuss Eq.Q that contains an inte- 
gration over the {K, u>) plane. This equation natu- 
rally displays a splitting of the heat flux into s- and 
p-polarized waves, and into propagating (K < lo/c) 
and evanescent waves (K > lo/c). The denominators 
account for multiple reflections through a Fabry- 
Perot term 1 — r z\ r z2 e ' 2lz d - The Planck function 
72 acts as a temperature-dependent frequency filter 
that cuts off frequencies much larger than ksT/h, 
i.e. beyond the near infrared at room temperature. 
As 73 ~ K for large K parallel wave vectors (deeply 
evanescent waves), there is also a wave vector fil- 
ter (e~ 273 <i ): wave vectors much larger than l/2<7 



do not contribute to the heat transfer at small gap 
sizes. This also implies that at sub-micron distances 
dcAj-, the evanescent contribution is much larger 
than the propagating one, leading to an enhanced 
heat flux. 

We show in Fig. 2 results obtained using a local 
dielectric constant. We consider a non-magnetic 
metallic medium characterized by a Drude model 
£ i,2(^) = £fc — ^p/i^ 2 + iuv) where e b accounts for 
the bound electron contribution, ui p is the plasma 
frequency and v is the damping coefficient. This 
model is appropriate for frequencies up to the in- 
frared range where the metallic response is mainly 
due to the conduction electrons. In this paper, we 
present results either for gold (e& = 1, u> v = 1.71 10 16 
S _1 , v = 4.05 10 13 s _1 ) or for aluminium (e b = 2, 
uj p = 2.24 10 16 s~\ v = 1.22 10 14 s - \ and we use in 
Sec. (3) vf = c/148 where vf is the Fermi velocity 
and c is the light velocity). 

Fig. 2 demonstrates that the increase of the heat flux 
levels off below distances of 10-30 nanometers, as was 
found in previous papers by Polder Q, Loomis [3] 
and Volokitin . The saturation is due to a strong 
s-polarized contribution. Only for distances below 
0.1 nm is the flux dominated by p-polarized waves, 
but in this regime, the local model is no longer valid 
(see Fig. 5 below). We note that in practice, with dis- 
tances in the nm-range, the s-polarized contribution 
dominates the heat flux. 

We now discuss the behaviour of the reflection coef- 
ficients in the (K, w)-plane (see Fig. 3). This points 
to the origin of the leading s-wave contribution. We 
plot the imaginary part of the reflection factors that 
is proportional to the heat flux [Eq.Q]. In partic- 
ular, also the local density of electromagnetic states 
(LDOS) is controlled by the imaginary part of the 
reflection amplitudes, as discussed in Refs. [HI [56}. 
First of all, we observe that lmr s (K,ui) covers a 
larger domain in the (K, w)-plane and takes larger 
values than its p-polarized counterpart. For the lat- 
ter reflection coefficient, one has at large K: 
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where R = (e b ~ e 3 ) / (e b + e 3 ) and u 2 sp = iv 2 /(e b + e 3 ). 
If medium 3 is vacuum and the background polariza- 
tion is negligible, R = 0, and the surface plasmon- 
polariton resonance occurs at oj sp = lo p /\/2. This 
resonance implies a peak in the near-field radiation 
spectrum (4|| as seen in Fig. 3. It lies for typical met- 
als in the UV, way above the frequency range that 
contributes significantly to the heat flux. Note that 



the asymptotics (0 becomes relevant only for ex? 
tremely large K- vectors where K ^> uo^Jei/c ^> uj/c, 
this is why the p-polarization becomes dominant 
only at very short distances (see Fig. 2). For the 
s-polarization, we have in the same range of K, 
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which tends to zero like 1/K 2 . This is the reason 
why the s-polarized contribution is often discarded 
when looking at the asymptotic behaviour [4l[ . But 
as shown on Fig. 3(a), there is a region where Im(r s ) 
has large values before decaying, corresponding to 
the wide interval uj/c <C K <g; ^/|ei| uj/c. We detail 
in the Appendix the behaviour of the reflection co- 
efficient and how to find the borders of the regions 
sketched in Fig. 3 (a). The result is an upper wave 
vector given by 



(7) 



Thus, we predict a saturation of the s-polarized heat 
transfer at gap distances smaller than 
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where the metal skin depth S is defined by l/5(ui) = 
(w/c)Im^/e. For frequencies between v and lo p , 
5 ~ c/lu p . For gold, the skin depth in this region is 
S = \^2c/u} p w 25 nm. It follows that the saturation 
distance is given by the skin depth at frequencies 
higher than v. We note that for gold, d m in — 18 nm. 
This is of the same order of magnitude as the cut- 
off distance in the experiment of Kittel et al. [27| . 
and, incidentally, also comparable to the electron 
mean free path. To summarize this section, we have 
found that the derivation of the heat flux between 
two metallic surfaces using a local dielectric constant 
predicts a saturation of the flux at a distance given 
by the skin depth. 



III. NEAR-FIELD RADIATIVE HEAT FLUX 
USING A NON-LOCAL MODEL 



We now turn to a non-local description of the heat 
transfer. There are several reasons to investigate the 
role of non-local effects in the heat transfer. First 
of all, non-local effects becomes significant at short 
distances. It has been shown that non-locality can 
explain the anomalous skin effect [47| | and has a very 
important effect on the lifetime of an excited atom 



or particle near a surface [4lj, [42], |43| . It has been 
seen that it has a si gnifi cant impact in the problem 
of near-field friction[48| . It has also been suggested 
that saturation of the heat flux could be due to non- 
local effects [27|. In addition, it is desirable to an- 
alyze the interplay between the skin depth found 
above and the mean free path. 

Temporal dispersion (i.e., frequency dependence of 
optical properties) appears when the electromag- 
netic (EM) field varies on a timescale comparable 
to the microscopic timescales of the medium where 
it propagates. A non-local behaviour (i.e. spatial 
dispersion or k-dependence of the optical proper- 
ties) is expected if the EM field varies appreciably 
on length scales given by the microscopic structure 
of the medium. 



ef M (k, u>)=e b - 2( \. {u, [f t (z, u) (10) 
w z (a> + iv) 

- Zz 2 h{z, u)]+iv [f t (z, 0) - 3z 2 fi(z, 0)] } 

where e b is the bulk contribution to the dielectric 
constant. It describes the interband contributions 
and it is constant in the following as these transitions 
do not play any role in the frequency range that 
we address. The Lindhard functions fij(z,u) have 
arguments z = fc/2fcp and u = (u> + iv)/kvp with 
Uf the Fermi wave vector, and are given by 

, . 1 1 - (z - uf z-u + 1 

fi{z,u) = - + In - 

2 8z z — u — 1 
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For metals, there are several microscopic length 
scales related to the Fermi velocity vf of the con- 
duction electrons. A first one is the electron mean 
free path vf/v, typically 20 nm for gold at ambi- 
ent temperature in the bulk. The second one is 
the charge screening length in a plasma of elec- 
trons called the Thomas-Fermi length, on the order 
of vf/^p- A third length is the Fermi wavelength 
1/kp = h/m*VF (where to* is the effective mass of 
the electron). It sets a lower limit for the spatial 
variations of the electron density in the metal and 
is often comparable to the Thomas-Fermi length. 
The fourth characteristic length is the distance vf/lo 
travelled by an electron during one period of an ap- 
plied EM field. This length governs an enhanced 
absorption by evanescent waves with K > uj/vf- 
This process is called Landau damping and consists 
in the creation of electron-hole pairs by absorption 
of photons. 

In order to account for the bulk effects we use 
two different dielectric functions: the Lindhard- 
Mermin (LM) and the Boltzmann-Mcrmin (BM) 
formulae [4l[. The LM dielectric function is also 
known, e.g., as the rando m p hase approximation 
(RPA) [4l|, Kliewer-Fuchs [49| constants or jellium 
ones. Other types of non-local dielectric functions 
are possible: the hydrodynamic model is an ap- 
proximation at small wave number [5fl | ; Feibelman's 
model {5lJ focusses on surface effects and has difficul- 
ties in taking bulk absorption into account, which is 
playing a significant role in heat transfer. We follow 
the notations of Ford and Weber for the longitudinal 
and transverse dielectric functions 41] 
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The limit u — > has to be taken with a positive 
imaginary part so that 

^ > 0) = i + i^ln|f±l| (13) 
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A semiclassical approximation of these formulae is 
obtained for wave vectors k much smaller than kp<, 
taking 2 = 0. This gives the Boltzmann-Mermin 
formulas 
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and 
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A few remarks are in order here. First, the Drude 
(9) formula is recovered at small k (u large and z small). 
Second, the variable u compares k to a combination 



of the mean free path vp/v and the distance covered 
by an electron during a period of the field Vf/oj, 
that can be considered as an "effective mean free 
path" (52j- Third, at very large wave vectors, the 
logarithms in Eqs. (|17ll8[) describe Landau damp- 
ing. Indeed, even for v = 0, they imply Im(e) > 
for k > u/vp 43]. Finally, it is seen that at very 
large wave vectors, there is a sharp cut-off in the 
imaginary parts of the Lindhard-Mermin dielectric 
functions: 
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Thus, fields oscillating with spatial periods smaller 
than half the Fermi wavelength cannot be screened 
by the electron plasma. 

We now account for microscopic surface effects that 
modify the reflection amplitudes. For the sake of 
simplicity, we use the infinite barrier model (also 
known as SCIB) that considers that electrons un- 
dergo specular reflection at the boundary [4l|. A 
model considering diffuse reflection of electrons is 
also available [53j. In our specular case, the reflec- 
tion coefficients are computed in terms of surface 
impedances as follows 



with 
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where under the integral, k 2 = K 2 +q 2 . K is the unit 
vector in the direction of the parallel wave vector K . 
As we account for spatial dispersion by using a non- 
local model, the reflection coefficients depend on id 



and if in a more complicated way than the Fresne^ 
formulas. 

One should note that in this approach, we do not 
tackle several effects that occur on the atomic (sub- 
nm) scale. The electron density, which is modified 
near the interface, is treated here with a step form 
and the addition of surface currents [ill, [4i| . Several 
authors [U [STj] showed that a self-consistent calcu- 
lation leads to a continuous variation of the elec- 
tron density between the bulk density and vacuum 
and that this can be described by an effective mean 
displacement of the surface, of the order of a few 
angstroms. Phenomena like electron tunneling also 
occur as the two surface approach each other on this 
scale and mutually influence their electron density 
profiles. We do not take this tunneling into account 
as it is clearly negligible in the nanometer range. 

We show on Fig. 4(a) the imaginary part of r p , 
at fixed u>. It is related to the local LDOS (see 
section IV). An interesting finding is that the lo- 
cal description leads to a plateau for large K (non- 
retarded approximation) that does not agree for 
any value of K with the non-local model. The lo- 
cal quasistatic approximation that has been often 
used, thus yields an incorrect value of Im(r p ) for 
a very broad range of frequencies. The curve la- 
belled 'longitudinal quasistatic' is based on neglect- 
ing the first term in Eq. (|25p . involving the trans- 
verse part of the dielectric function. We see that 
this term nevertheless contributes at wave vectors 
K < 1/5(uj) w 10 s m _1 . For larger K, the non-local 
calculation leads to an increase of Im r p by roughly 
one order of magnitude, that we attribute to Landau 
damping. Finally, we observe that for wave vectors 
larger than k F ~ 10 10 m _1 , the non-local models 
predict a strong decay of Im(r p ) as compared to the 
local model. 

The s-polarized reflection coefficient Im(r s ) is plot- 
ted in Figure 4(b). Differences to the local calcula- 
tion are barely visible in the domain K < 5 10 s m _1 
where Im(r s ) takes significant values and contributes 
to the heat transfer. We thus expect only small cor- 
rections to heat transfer from the non-local models. 

Figure 5 presents the heat flux as a function of the 
gap distance. We display the fluxes due to s and p 
polarizations when using both a local and the two 
non-local models introduced above. Although the 
validity of the models is questionable for distances 
smaller than 1 nm, we display the flux at smaller 
distances in order to analyze their physical content 
when d — > 0. What is important here is that the 



local and non-local heat fluxes arc identical up to 
distances on the order of the Thomas-Fermi length 
vp/ujp. ft appears that the small modifications of 
Im(r s ) give the same final result after integration 
over K and u>. A small increase of the heat flux 
due to the onset of Landau damping is observed in 
the p-polarizcd contribution, but in a regime where 
s- waves dominate and level off. Another observation 
is that the two non-local models are superimposed, 
showing that the Thomas-Fermi length is sufficient 
to describe the large K decay of the dielectric con- 
stant. Finally, at very short distances (below the 
Thomas- Fermi length or the Fermi wavelength), the 
non-local models remove the 1/d 2 regime of the p- 
polarized flux. 

We now illustrate how the non-local models sup- 
press this 1/d 2 dependence. We have plotted in Fig. 6 
the p-polarized contribution to the heat flux in the 
(K, aj)-plane, but removing the decay term e -21 ™' 73 -' 
and the Planck function /°(T) that act as filters. 
What we plot is thus lm{r^) 2 /\1 - {rf) 2 e~ 2 ^ d \ 2 . 
Figs. 6 show a locus that follows the dispersion re- 
lation of the surface plasmon-polariton. ft is seen 
that it has two branches [45|, [54|. They split at 
a wave vector of order 1/d that is pushed towards 
large K as the gap size is decreased. When non- 
locality is included, the flat asymptote at frequency 
uj sp = ujp/\/2 for large values of K becomes disper- 
sive and approaches lu = vpK in Fig. 6(b). But what 
is important here is that the far IR branche of the 
resonance will not be able to be shifted to the large 
K region when the gap size decreases because of the 
cut-off at w = vfK. This removes the divergence 
of the heat flux due to the p-polarized evanescent 
contribution in Eq|3]when d — > 0. ft provides an in- 
trinsic cutoff at large K which is different from the 
distance d. 

The main conclusion of this section is that the lo- 
cal calculation is in practice sufficient when comput- 
ing heat fluxes between two metallic surfaces a few 
nanometers apart. The second conclusion is that 
non-locality removes the universal heat flux diver- 
gence at short distance as expected. 



IV. DISCUSSION AND CONCLUDING 
REMARKS 



In this last section, we try to gain some insight on 
the physical mechanisms responsible for the near- 
field heat transfer in s-polarization between two par- 



allel interfaces. We have plotted in Figure 7 th§ 
LDOS (irl |55| near a metallic/vacuum interface in 
vacuum. We recall that the local density of energy 
is the product of the LDOS by the mean energy of 
an oscillator given by hu>/(e k B T — 1). The LDOS 
is split in four contributions: magnetic and electric 
fields, s- and p-polarization. For instance, the con- 
tribution of the evanescent s-polarized magnetic field 
to the LDOS is given by 

pf (z, u ) = Pv / + °° *K*Lf(K,u)Tm(r.) 
Ju/c 2|7 3 | u 

(25) 

where p v (uj) = lu 2 /it 2 c 3 is the vacuum density of 
states and f(K,uS) — 2(^j-) 2 — 1. Again, the prop- 
erties of the material control the LDOS via Im(r s ). 
Fig. 7 shows that the propagating terms are neg- 
ligible. Furthermore, the leading contribution in 
the infrared (u) ~ I0 13 . . . I0 15 s _1 , where the room- 
temperature thermal spectrum peaks) is clearly due 
to s-polarized magnetic fields. It follows that a 
metallic half-space generates a very large magnetic 
energy in a vacuum close to the surface. This quan- 
tity is relevant to analyze the heat transfer through 
an interface. Indeed, as the magnetic field is con- 
tinuous through an interface with a non-magnetic 
material, the magnetic field penetrates without re- 
flection. 

The large value of the magnetic density of en- 
ergy due to s-polarized waves near a metallic 
interface has been discussed recently [H, |5(|. 
Whereas the ratio c|i?|/|i?| takes a fixed value of 
1 for propagating waves, it becomes frequency- 
dependent for evanescent waves (K/k > I). For s- 
polarized evanescent waves, one can show using the 
Maxwell-Faraday equation that this ratio is given 
by yj f(K, lj) ~ \/2K/ko- Magnetic fields dominate 
in s-polarization. For p-polarized waves, the oppo- 
site trend |S|/c|B| ~ yj f{K, uj) is found, showing 
that electric fields dominate. If we want to know 
which of the magnetic s-polarized waves or the elec- 
tric p-polarized waves give the leading contribution 
to the LDOS, we have to compare the products 
f(K,uj)lm(r s ) and f (K , u)Im(r p ) . As we have seen, 
the s-polarized reflection coefficient is larger than 
than Im(r p ) for a metal at infrared frequencies and 
below, so that finally, the LDOS is dominated by its 
s-polarized magnetic component as seen on Fig. 7. 

It follows that retardation plays a key role as ob- 
served in Ref. (H)]. Accordingly, the heat trans- 
fer between a metallic nanoparticle and a half space 
(l6l . HH, H3| must be revisited accounting for mag- 
netic energy. It will be shown that the magnetic 



dipole yields the leading contribution [56(. 

The large magnetic fields can be traced back to the 
current density in the material. In s-polarization, 
the electric field E is tangential to the metallic in- 
terface and therefore continuous. It drives a sur- 
face current flowing within the skin depth 8, with 
an amplitude roughly given by aE. This suggests 
the following mechanism for the heat transfer be- 
tween metallic surfaces: fluctuating currents flow- 
ing parallel to the interface within the skin depth in 
medium 1 generate large magnetic fields at IR fre- 
quencies. These fields penetrate into medium 2 and 
generate large eddy currents which are dissipated by 
the Joule effect. In other words, radiative heat trans- 
fer in the near field is similar to nanoscale induction 
heating at infrared frequencies. 

We have seen in section 2 that the skin depth plays 
a key role (57j . The above argument provides a sim- 
ple picture for the phenomenon. The skin depth 
depends on the frequency. We stress that the cutoff 
distance seen by Kittel [27J and that we found above 
is linked to the skin depth evaluated at the frequen- 
cies contributing to the largest parallel wave vectors, 
Lomi. For gold, this skin depth is 8 = \Z2c/uj p 
25 nm. Our analysis leads to a number of predictions 
that should be measurable. Measurements of the 
heat transfer such as reported by Kittel should be 
able to detect the skin depth dependence by chang- 
ing the metals. As seen on Figure 8, the plasma 
frequencies of a number of metals are not very dif- 
ferent. They all give (local) cut-off distances in the 
range of 10 to 200 nm. The differences should be 
measurable. A material like cobalt is expected to 
saturate at larger distances than metals like cop- 
per, gold or aluminium. Interestingly, cobalt could 
also be a test-case study for the saturation due to 
non-locality as the p-polarized contribution becomes 
larger than the s-polarized contribution near 1 nm. 
Another interesting issue is the heat flux between 
two different metals. We expect a saturation dis- 
tance governed by the smallest skin depth due to the 
product Im(rj 1 )Im(rg 2 ) in the heat flux formula. 

To summarize, we have shown that the radiative 
heat flux between two parallel metallic surfaces sat- 
urates when the gap size reaches a distance equal to 
the skin depth at a frequency equal to v. We have 
shown that the leading contribution to the flux is 
due to eddy currents generated in the medium. The 
non-local effects have been studied. They do not sig- 
nificantly affect the s-polarized fields but introduce a 
cut-off in the K dependence of the p-polarized fields. 
This cut-off removes the 1 /d 2 dependence of the flux 



at short distances. As the s-polarized fields domi- 
nate the heat transfer between metallic surfaces, the 
non-local corrections are negligible. Finally, we ob- 
served that the cut-off distances seem to be in the 
range of 10 to 200 nm for many metals. 



APPENDIX 



We explain in this section how we estimate the lim- 
its of the domain in the (K,u>) plane where Im(r s ) 
contributes to the heat flux. As is shown on Figure 
3, the (K, u>) plane can be divided into four areas. 
Point A is the intersection of the 4 borders. In all the 
cases, we consider only evanescent waves: K 3> &o 
with ko = uj/c. 

We address first the division of the (K, to) plane be- 
tween large K and smaller values. This underlines 
the different behaviours of regions 1 and 3 on one 
hand, and regions 2 and 4 on the other hand. The 
perpendicular wave vector 71 is given by 

K 2 + 1 2 1 = £l fc 2 (26) 

where fco — w/c. This shows that we have two 
regimes. To leading order, we have 7 2 ~ —K 2 at 
very large K (regions 2 and 4) and 7 2 ~ e\k 2 at 
smaller K (region 1 and 3). The transition occurs 
at a critical wave vector K 2 ~ | e 1 Ajq | . This gives a 
critical wave vector given by 

^H^V^Mfco^y^T ( 27 ) 

where the last equality applies to the Drude model 
at frequencies u> <C u> p /yfet>. Values of r s in both 
regimes are now given. To leading order, one finds 

f -1-2-L— - (region 1,3) 
r<« , 2 ^ lko (28) 

I 4iP (£l " 1} (regi ° n 2 ' 4) 

At large K, Im(r s ) decreases to small values that do 
not contribute significantly to the heat flux integral. 

We address now the horizontal division of Fig. 3. The 
upper region is given by domains 3 and 4 and the 
lower one by domains 1 and 2. This limit is due to 
the different behaviours of e(u>) if lo <C v (domains 
1 and 2) or w > 1/ (domains 3 and 4). The first two 



asymptotic orders are 



iul/wv £ (region 1,2) 



(29) 



(region 3,4) 



The low-frequency expression is also known as the 
Hagen- Rubens formula. In Table HI we give the cor- 
responding asymptotics for Imr 5 in the four regions. 



Region 


Characteristics 


Im(r s ) 


1 


far IR, small K 


V2vc K 
Up \/u 


2 


far IR, large K 


2 

p u 
4,vc 2 I? 1 


3 


near IR, small K 


vc K_ 

Up UJ 


4 


near IR, large K 


"1" 1 a 



TABLE I: Asymptotic behaviour of Im(r s ). A local 
Drude model is taken for e(u) with plasma frequency 



L0 P and relaxation rate v. 



As a function of frequency, the critical wave vector 



behaves like K c w (lj p / c)(lj / v) 1 / 2 in the far infrareS 
(small frequencies) and like K c w lo p /c for larger 
frequencies. These two lines cross at u> ~ v which 
is the point A marked in Fig. 3. At this point, the 
imaginary part of r s {K, oj) reaches its maximum. 

According to Eq. ([2T|) . Im(r s ) takes significant values 
for K lower than K c — lo p /c. This limit yields a sat- 
uration length 1/K C = c/oj p . Note that this length 
is related to the skin depth as S = - — , } , , ~ —t= — . 

Im(VeiKo) V2uip 

At low frequencies (regions 1 and 2), e is purely 
imaginary, leading to S ~ \^2/K c , while in the high 
frequency regions 3 and 4, S ~ l/K c . Hence, at each 
frequency, the cutoff wave vector is essentially given 
by the inverse skin depth. 
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